SPECIFICATION 
TITLE 

"METHOD FOR AUTOMATICALLY SETTING AN X-RAY DOSAGE FOR 
PRODUCING AN X-RAY TOMOGRAPHIC IMAGE" 

BACKGROUND OF THE INVENTION 

Field of the Invention 

The present invention is directed to a method for automatically setting an X- 
ray dosage for producing a tomographic image, using an X-ray CT device. 
Description of the Prior Art 

In the production of images of an examination subject in a computed 
tomography system, it is unavoidable that a certain level of noise will be present in 
the image. In general, the level of noise can be reduced by increasing the tube 
current, and thereby increasing the X-ray dose, however, obviously this cannot be 
done indiscriminately because overexposure of the examination subject to X-rays 
must be avoided. 

Moreover, different portions of a human examination subject exhibit different 
attenuations, due to respectively different compositions and physical dimensions. 
Particularly in computed tomography examinations wherein a larger portion of the 
body is scanned, which may encompass such differently dimensioned body regions, 
it is undesirable, from the point of view of image quality, to select one and only one 
tube current/X-ray dose for the entire scan, because the dose may be suitable for 
scanning one portion of the subject, but may not be suitable for scanning another 
portion. If a single dose is employed for the entire scan, it must be selected at a 
compromise value. 

Various investigations have been undertaken to analyze the relationship 
between the level of noise in a CT image, and the respective values of various 



parameters that are associated with the production of the CT image. One example 
of such an investigation is the article "Noise Due To Photon Counting Statistics In 
Computed X-Ray Tomography," Chesler et al., Journal of Computed Assisted 
Tomography, Volume 1, No. 1 (1977), pages 64-74. 

SUMMARY OF THE INVENTION 
It is an object of the present invention to provide a method for adjusting the X- 
ray dose in the production of a computed tomography image wherein an overall 
image with constant image quality is produced without overexposing the examination 
subject to X-rays. 

The above object is achieved in accordance with the principles of the present 
invention in a method wherein the tube current is limited automatically at higher 
levels of attenuation of the examination subject, while still insuring appropriate image 
quality. A tube load computer, preferably mounted on the gantry of the CT 
apparatus, employs values obtained in reference measurements for different 
attenuation values and, dependent on the measured attenuation and the particular 
organ or body region being examined, the tube current is appropriately limited. The 
reference measurements can be obtained using a water phantom. On the basis of 
water phantoms having different attenuation, the interrelation of the desired image 
noise, tube current, maximum attenuation of the subject, and detector voltage are 
determined. These values are stored in a table for a particular tube voltage and slice 
thickness. The table is consulted by the load computer and operating parameters 
are set for examination of a particular organ or body region by selecting, from the 
table, the parameters for the water phantom having an attenuation most closely 
approximating the attenuation of the organ or region under examination. 
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DESCRIPTION OF THE DRAWINGS 

Figure 1 is a schematic illustration of a CT apparatus for use in explaining the 
inventive method. 

Figure 2 is a flowchart illustrating the sources of noise in a CT image. 

Figure 3 is a graph showing the detector voltage for various modes of 
operation of a CT apparatus. 

Figure 4 is a graph showing the ratio of the actual tube current to a reference 
tube current as a function of attenuation. 

Figure 5 is a graph showing noise as a function of the subject attenuation. 

Figure 6 shows the attenuation-dependent tube current for use in the 
inventive method. 

Figure 7 schematically illustrates adjustment of the tube current in accordance 
with the invention for scanning different body regions. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

First the theoretical interrelation between measurement system, absorption 
and noise is summarized and is applied to measurements on the Plus 4. 
Subsequently, procedures for automatic dosage are discussed. 

The number of quanta NO emitted by the tube during the measurement of a 
projection and of the slice fade-in h depends on the tube high voltage kV, the tube 
current IRohre, the pre-filtering V and the wedge filter W: 

N 0 =N 0 (kV, l R 6hre, h, V, W) (1) 

For a system as shown in Figure 1 , the pre-filtering can be assumed to be 
constant. The wedge filter causes a detector-channel-dependent variation of the 
quanta N 0 . The high voltage, as well as the pre-filtering and wedge filter, determines 
the spectral energy distribution of the x-ray spectrum. 
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For the functional interrelation of N with the high voltage U, the tube voltage I, 
and slice thickness h, the following holds: 



N = N 0 * 



( jj A 2 - 26 r r f \ 



h_ 

[to ) [ h o 



(2) 

It can be seen that N ~ I, N ~ h and N ~ U 226 . If U = const, h = const, then the 
number of quanta depends in a linear fashion on the tube current I. 

After passing through a homogenous medium having thickness d and 
absorption jli, the following holds for the number of quanta N: 
N - No * e ^ d (3) 
The subject attenuation A is given as 

A=N 0 /N (4) 
N is the number of quanta that impinge on the detector and produce the signal U D et 
in a channel of the detector: 

U Det = C * N (5) 
The constant c is a property of the detector that describes the absorption of 
the quanta and their conversion into an electrical signal (Signal/Quant). C depends 
on the spectral energy distribution of the x-ray spectrum and on the spectral 
sensitivity of the detector. Without an absorber, N 0 is to be used instead of N in 
equation (5). 

The detector signal from equation (5) varies due to the quantum noise. In 
addition, the electronic noise U E must be taken into account. Thus, the following 
holds for the detector signal: 

Uoet= UQuanten + Ue (6) 
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The quantum noise is based on the Poisson statistic, so that the following 

holds: 



In connection with equation (4), the following holds for the quantum noise U Q after 
passage through a subject having attenuation A: 



According to equation (5), U 0 is the detector signal without absorber. From equation 
(8), it can be seen that there is an attenuation A 0 at which the quantum noise is equal 
to the electronic noise: 

Ao = (Uo/Ue) 2 /No = (c/Ue) 2 * No (9) 
Equation (6) can then be transformed to 

U D e, = c*No/A + c*jNo/A -KJe (10) 
The detector signal thus depends on system-specific quantities (pre-filtering, 
wedge filter, detector constant c) that are equal for all modes of operation of the 
scanner, as well as on quantities (kV, Ipahre, h) that are selected by the user, and of 
course also on the attenuation A of the subject of measurement. In addition to the 
noise of the measurement signal, which enters into the reconstruction via the line 
integral, the noise from the reconstruction itself must be taken into account as shown 
in Figure 2. 

For linear interpolation, the image noise a of a homogenous water disk in the 
center of the measurement field is calculated as follows: 



Uguanten = C * (N ± V# ) = 



c*N 



(7) 




(8) 



500 



(11) 



Mja^NEQ 
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The relative noise factor q r takes into account the noise of the reconstruction 
algorithms. According to the definition, q r = 1 for the Shepp-Logan core kernel with 
linear interpolation. a D is the scanning raster determined by the detector and related 
to the center of rotation. For the noise equivalent quanta, NEQ, the following holds: 



N * N 
NEQ= R 0 



a D * A 



(12) 



N R is the number of readings per scan (number of projections * reading number, e.g. 
1.5 sec: 1056 projections * 2). From equations (11) and (12), the relation can be 
seen between the number of quanta (dosage), subject attenuation and image noise. 

A 



500 



a * N 



(13) 



The following table shows the results of a measurement of the 20 cm water 
phantom in the Plus4 (UFC detector) mode (Quality/Noise): 



Parameter 


Value 


Derived quantity 


Mode 


RED 




KV 


140 




mA 


146 




Tl 


1.5 


Projections 1056 Reading number N R = 2 


SL 


10 




Kernel 


SP90 


q r =1 


FOV 


250 0/0 




Result 


cri40kv = 8.76 





With equation (13), N 0 can thus be calculated (a D =0.34 mm) : N 0 = 1.68-10 6 
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For determining the electronic noise: 



Measurement 


Evaluation 


Parameter Value 

Mode STATIC 
KV 140 
MA 146 
Tl 1 .5 sec 
SL 10 mm 
SP90 
FOV 

Storage FPA 

Control platform X-ray OFF 
Anode Red OFF 

Result 


Analyze 

Load 00001. RCT 
Row/Proj. Step: 2 
Format uV 
STAT / SD 
WA (19) 

Column/Channel lb = 350 
ub = 420 
Row Stat ~>MEAN 

Average electronic noise Ue 
in uV 



From the above values, U E = 32.93 uV 

For determining the strength of the detector signal: 



Measurement 


Evaluation 


Parameter Value 
Absorber WITHOUT 
Mode RED 
KV 140 
mA 146 
Tl 1 .5 sec 
SL 10 mm 
SP90 
FOV 

Storage FPA 

Control platform X-ray ON 
Anode Red ON 

Result 


Analyze 

Load 00001. RCT 
Row/Proj. Step: 2 
Format uV 
STAT/Mean 
WA(18) 

Column/Channel lb =350 
ub = 420 
Row /STAT — >MEAN 

j Signal strength U 0 in uV 



From these values, U 0 = 4460736 uV 



From equation (8), using the measurements it is thus possible to determine 
the quantum noise, and also the detector constant c, for each attenuation if a 12000- 
fold attenuation is assumed, U Q = 41 u.V results. 
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From equation (5), using the measured U 0 (sec. 4.3) and the value N 0 
determined in sec. 4.1 , the constant c can be determined as c = 2.65 jxV/Quant 

According to equation (10), Uoe: = c*N/A holds for the detector signal 
without noise portions. N can be calculated according to equation (2). The rotation 
time Tl can unproblemmatically be added to the equation. In 4.5, the constant c has 
been determined for the Plus 4 UFC detector. In Figure 3, the strength of the 
detector signal is shown for various subject attenuations (water phantom) and 
measurement protocols. It is clear that for 1 mm slices, beginning with a subject 
attenuation of 4000 (corresponding to a water phantom 40 cm in diameter), the 
electronic noise is greater than the detector signal. The attenuation in the shoulder 
can reach a value of up to 10,000. Figure 3 shows the detector voltage for various 
modes. 

As can be seen from Figure 2 noted earlier, the noise is composed of two 
portions: the noise from the measurement system (quanta and electronic noise) and 
the noise from the image reconstruction. The latter depends on subject properties 
(A, ^ w ), system parameters (ao), recording parameters (N R , N 0 ) and parameters of 
the image reconstruction (interpolation, kernel [q r ]). Subject properties and system 
parameters cannot be influenced; the parameters of the image reconstruction are 
normally fixed by the organ mode used, so that N 0 , i.e. the tube current, remains as 
a variable quantity that has to be optimized online during the measurement. The 
matching of the tube current to the maximum subject attenuation during rotation is 
necessary. 

The determination of the subject attenuation is already part of the dosage 
modulation and can be taken over. For a single line (or row) detector, this involves: 
1 . Offset correction of the raw data 
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2. Smoothing of the projections (e.g. sliding [or: variabl , flexible] averaging) 

3. Norming of the data (U M onitor/U k ), U k = Signal of the k-th detector channel 

4. Calculation of the maximum of all detector channels. 

The optimal dosage is determined by means of reference measurements. 
Since Equation (11) was calculated for round, homogenous water phantoms, at first 
only such phantoms need be considered. On the basis of water phantoms having 
different attenuation, the interrelation of the desired image noise a, tube current lR e f, 
maximum attenuation A of the subject and detector voltage UDet are determined. 
These values, U De t(A, a, I), are stored in a table for a particular tube voltage U k v and 
slice thickness h. 

If a is fixedly predetermined, i.e., the kernel is not modified, then it is sufficient 
to consider the detector signal U De t. From equations (6) and (7), the following 
results: 

No 

Uoet = C * + UE (14) 

A 

For an acceptable image quality, the electronic noise must be negligible. Dependent 
on the attenuation A, a minimum tube current I is to be predetermined, so that c*No/A 
> Ue holds. Udet is then simplified to 

No 

UDet = C* (15) 

A 

Ul e J ,a ,N* e/ ,A Re/ are the parameters for the reference measurement of a phantom 
having attenuation A, so that for a given dosage N 0 an image noise a and a detector 
signal U De t is produced. If a subject having a different attenuation A act , i.e., a 
different diameter, is now measured, then N 0 act is to be adapted so that the following 
holds: 
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T j Re f , a j j act 

U Det ~" U Del 



(15) 



Applying equations (2) and (15), the following results: 



*Re/ 



2.26 



jjkV 



h 



(16) 



In addition, equation (16) takes into account the tube voltage and the slice thickness. 
If these are not modified, the new tube current l act can be calculated: 



/ — «£<_ * / 

act . L Re f 

A Re/ 



(17) 



Figure 4 a was calculated corresponding to equation (17), whereby a round 
water phantom having a diameter of 20 cm (A = 46.5) was used as a reference, and 
norming took place to l re f = 146 mA (see sec. 4). It is clear that for phantoms having 
a diameter of greater than 25 cm, corresponding to an attenuation of 122, the 
required tube current can no longer be set! In the shoulder region, the attenuation is 
greater than 10000! It is therefore necessary to permit a higher degree of noise as 
the attenuation increases. In Figure 5, the noise is plotted over the attenuation, if the 
tube current is held constant (blue curve) and the tube current according to Figure 4 
is increased. 

Because the maximum tube current l max is predetermined by the system, it is 
necessary to ensure that 7 f continues to hold. However, in clinical 

act max 

application the outlook is not as negative as it appears according to Figures 4 and 5. 
It will not be expected to examine a shoulder with the low degree of noise of an 
abdominal measurement. 

The tube current can be limited by reducing the gradient of the change of the 
momentary tube current with increasing attenuation. In Figure 6, this is achieved by 
dividing the attenuation region into three segments having decreasing gradients. 
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Comparison with Figure 4 shows that the increase in the tube current is reduced 
from a factor of 120 to a factor of 5. The division of the attenuation can ensue in a 
further region having respectively different gradients, so that there results from 
equation (17): 

h« = giorgan, A organ * * / Re/ (18) 

A Re/ 

g (organ, A org an) is predetermined dependent on the organ mode and the maximum 
attenuation, if this is known for example from topogram measurements. 

F user is a factor permitting the user to adjust the reference current l Re f to his 
preference in order to achieve a lower degree of noise or a higher savings in dosage. 
The factor should be adjusted once and then left unchanged. 

As in dosage modulation, in automatic dosage the tube current is also 
modified online during the measurement. In contrast to modulation, which in the 
present implementation only reduces the tube current, in automatic dosage it is 
necessary to increase the current. 

One approach is shown in Figure 7, wherein the automatic control in 
accordance with the invention is designated "auto." As was also previously the case, 
the tube load computer takes into account the nominal currents required by the user 
for the selected organ program, and interactively carries out an optimization process 
in order to avoid cooling pauses. If the maximum thermal tube load is not achieved 
using the prescribed scanning protocol, the load computer simulates possible 
maximum limits. In this way, the nominal currents of all scans can be increased 
proportionally, or the magnitude of the increase of the currents takes place in an 
organ-specific manner. These limits are downloaded for control purposes, and 
predetermine the maximum range of variation for the tube current. 
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The following table describes the scans shown in Figure 7: 



Curr nts 


Remark 


Scan 1 


Scan 2 


Scan 3 


'nom 


Selected by user 


I i 

■nom 


1 2 

mom 


1 3 

mom 


'tLC max 


Calculated by 
TLC 


ItLC max 


ItLC max 


ItLC max 


I auto 


Tube current of 
the automatic 
dosage 


lauto 1 (Rotation) 


lauto 2 (Rotation) 


lauto 3 (Rotation) 



Table 1 



lauto 1 (Rotation)<= l-rLCmax' i = scan number (18) 

An online tube load computer can be located on the gantry, which monitors 
the tube load during the measurement and limits the momentary tube current, so that 
all planned scans can be measured with a "reasonable" tube current. It is to be 
avoided that, for example, two of three spirals are measured with an optimal tube 
current, while the current must be very strongly reduced for the third. 

In examinations using contrast agents, a scan region is measured twice. The 
control unit notes the maximum attenuation and the modulation function and applies 
them in the second scan. If necessary, interpolation can take place in the z 
direction. 

The maximum attenuation is determined ahead of time from two orthogonal 
topograms. Using the algorithm that is also used in the control unit for the 
interrelation between attenuation and tube current, a good estimation of the required 
tube current can be produced ahead of time. 

Although modifications and changes may be suggested by those skilled in the 
art, it is the intention of the inventors to embody within the patent warranted hereon 
all changes and modifications as reasonably and properly come within the scope of 
their contribution to the art. 
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